Bose-Einstein condensation of trapped interacting spin-1 atoms 
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We investigate Bose-Einstein condensation of trapped spin-1 atoms with ferromagnetic or an- 
tiferromagnetic two-body contact interactions. We adopt the mean field theory and develop a 
Hartree-Fock-Popov type approximation in terms of a semiclassical two-fluid model. For antifer- 
romagnetic interactions, our study reveals double condensations as atoms in the \m,F = 0) state 
never seem to condense under the constraints of both the conservation of total atom number N and 
magnetization M. For ferromagnetic interactions, however, triple condensations can occur. Our 
results can be conveniently understood in terms of the interplay of three factors: (anti) ferromag- 
netic atom-atom interactions, M conservation, and the miscibilities between and among different 
condensed components. 
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I. INTRODUCTION 

Research in dilute atomic quantum gases remains one 
of the most vibrant areas in physics almost ten years after 
the landmark discovery of atomic Bose-Einstein conden- 
sation (BEC). Increasingly, new experiments are reveal- 
ing the rich possibilities afforded by the internal elec- 
tronic state structures of an atom, e.g., the generation 
of atom entanglement in the Mott phase through con- 
trolled collisions from the relative displacement of the 
optical lattice potentials of the respective atomic inter- 
nal states 0] and the recent observation of condensation 
of fermionic atom pairs |2(. 

It has been known for a long time that the spinor de- 
grees of freedom of an atom becomes accessible if a far- 
off-resonant optical trap is used to provide equal confine- 
ment for all Zeeman states, instead of the more widely 
used magnetic traps for spin polarized atoms 0, 
Several earlier experiments have produced fascinating ob- 
servations of spinor condensates, a superfluid with inter- 
nal degrees of freedom, e.g., with 23 Na atoms in F = 1 
and F = 2 and 87 Rb atoms in F = 1 H 03 and 
spin domains and interdomain tunnelling 

S lhi . as well as the generation of coreless vortex states 
llq . IrH ITsf . These properties exist only because of 
the spinor nature of the condensate order parameter, and 
thus are generally not expected to occur in a magnetically 
trapped condensate. 

Despite these and other related successes with spinor 
condensates, our knowledge remains limited regarding 
the condensation thermodynamics of spin-1 atoms. In a 
sense, the spin-1 condensate constitutes a type of quan- 
tum fluid unfamiliar to many of us. On the experimental 
side, it remains a significant challenge to produce a spinor 
condensate, as evidenced by the disproportionally small 
numbers of spinor BEC experiments in operation. In this 
article, we reconsider the topic of the condensation ther- 
modynamics for a system of trapped spin-1 atoms. Of 



particular interest to us is the question of the so-called 
double condensations for a spin-1 system constrained by 
two global conservations [jjj- Using the Bogoliubov- 
Popov approximation, Isoshima et al. first investigated 
the thermodynamics of the BEC phase transition for a 
spin-1 gas [l||. Huang et al. studied analytically the ef- 
fect of a magnetic field on the transition temperature |20| . 
While an attempt to find the zero-magnetic-field phase 
diag ram was made through numerical simulations in Ref. 
[ijj , there still exist several question marks to the overall 
picture of BEC for a spin-1 Bose gas, especially for ferro- 
magnetically interacting atoms such as 87 Rb. Limited by 
the computation procedure within the Bogoliubov-Popov 
approximation, only a few data points were made avail- 
able in the earlier studies by Isoshima et al. 19] . The 
lack of focused experimental efforts also indirectly dis- 
couraged a detailed investigation of this problem until 
now. 

In this paper we systematically investigate the phase 
diagram of a spin-1 Bose gas for both ferromagnetic and 
antifcrromagnctic interactions. Instead of the Bogoli- 
ubov approach, we will adopt the Hartree-Fock-Popov 
approximation and employ a semiclassical approximation 
to the noncondensed atoms within the mean field the- 
ory. We will also enforce the thermodynamics for a finite 
trapped system with a fixed total atom number N and a 
total magnetization M. Recent studies have significantly 
verified the accuracy of this approximation when applied 
to similar systems |2J . As we will illustrate in this work 
our results indicate that double condensations will oc- 
cur for a spin-1 gas with antiferromagnetic interactions, 
while triple condensations are more likely for ferromag- 
netic interactions. 

This paper is organized as follows. In Sec. II, we re- 
view the additional features of a BEC for an ideal gas 
of spin-1 atoms. This is followed by the discussion of 
an interacting spin-1 gas in Sec. Ill and a brief sketch 
of the Hartree-Fock-Popov theory used for our investiga- 
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tion. We outline the detailed numerical algorithm used 
to solve the coupled two-fluid model quantum gas at dif- 
ferent temperatures in Sec. IV and present the results of 
our study in Sec. V. We conclude with some discussions 
and remarks in Sec. VI. 



the temperature is lowered, the |+) component first con- 
denses for a system with a positive M because its phase 
space density is largest, reflecting the fact that N + is the 
largest component population. Thus we first arrive at 
/i + = 0. This consideration leads to the critical temper- 
ature T\ governed by the following equations: 



II. BEC OF AN IDEAL GAS OF SPIN-1 ATOMS 

In this section, we briefly review the phenomenon of 
a BEC for a trapped noninteracting gas of spin-1 atoms 
following the pioneering study of Isoshima et al. |19| . 
At thermal equilibrium, we adopt the standard Bosc- 
Einstein distribution, and treat the spinor degree of free- 
dom as degenerate internal states in the absence of an 
external magnetic field. The average number of atoms at 
each single atom state of an energy £j for the component 
\F = 1, mjr = i), i = +1,0, and —1 (hereafter is 
then conveniently given by 
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with f3 = 1/(&bT) at temperature T. ks is the Boltz- 
mann constant. The fugacity can be expressed in terms 
of the chemical potential for the ith component [ii as 
i-i = exp(/3/x;). In the thermodynamic limit, one can fol- 
low the usual approach by making a semiclassical approx- 
imation for a continuous description of the single particle 
density of states, and treating the ground state popu- 
lation separately as it can become macroscopic due to 
Bose-Einstein condensation. The total number of atoms 
for a given internal state in all motional excited states of 
the trap is thus found to be 



N = 



M 
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Ml)+53(e^)+ff 3 (e 2 ^)], (3) 
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On further lowering of the temperature, however, the re- 
maining two components |0) and |— ) condense simulta- 
neously, rather than sequentially with the less populated 
component of the two condensing last. This is precisely 
due to the conservation identities as discussed before. 
The relationships fi± = fj, ± r\ and fiQ = fi lead to a 
mathematical certainty: when fi + is zero, if either /j,q 
or becomes zero, both must be zero. At this sec- 
ond critical temperature X2, both (1q = and /i_ = 0, 
which imply that the |0) and |— ) components condense 
simultaneously This second condensation where all three 
components condense, occurs at the temperature T2 of 
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In Fig. ^ we illustrate the M dependence of the dou- 
ble condensations for an ideal Bose gas of spin-1 atoms. 
T c = [N/g 3 (l)] 1 / 3 (hLj/kB) « OMN^^uj/kB is the con- 
densation temperature for M = N, i.e., for a single com- 
ponent gas with all atoms polarized in state |+). 
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where we have assumed atoms are confined in a spheri- 
cal harmonic trap with a frequency u> independent of the 
atomic internal state \i ). gtjx ) = Y^=i( xn / ' n ^) * s the 
standard Bose function [22,|23|]. We note that the conser- 
vations of the total number of atoms N = N + + Nq + N- 
and total magnetization M = N + — iV_ lead to the chem- 
ical potentials for different spin components expressible 
as n± = fx ± 77 and /j,q = [i. These identities remain 
valid in the presence of atom- atom interactions, /i and 
•q are effectively independent Lagrange multipliers used 
to guarantee the conservation of N and M, respectively. 
Taking the single atom trapped ground state to be zero 
energy, the Bose distribution (JIJ shows that is nega- 
tive at high temperatures and reaches zero when the spin 
component |i) condenses. r\ is positive (negative) for a 
positive M (negative), which acts as a fictitious applied 
magnetic field physically. 

As was first pointed out in Ref. [n|, there exists an 
interesting double condensation phenomenon for a spin- 
1 gas because of the presence of M conservation. When 



III. BEC OF AN INTERACTING GAS OF 
SPIN-1 ATOMS 

A. Formulation 

Our model system of the interacting spin-1 atoms is 
described by the following Hamiltonian in second quan- 
tized form 



H 



V 2 + V ext (r) 



+ 



°j4 (F a ) l3 ^l (F a ) kl ^} 



(6) 



where ipi(r) [ipl(r)] is the quantum field for annihilating 
an atom in state |i) at r, and i,j, k, I = +, 0, — and a = 
x, y, z. Repeated indices are assumed to be summed 0, 
Q. F a=Xt y tZ are spin-1 matrices given by 
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FIG. 1: BEC for a gas of spin-1 atoms with M > (M < 0). 
For noninteracting atoms, the |+) (|— )) component condenses 
first at T\ (dashed line) while the |0) and | — ) (|+)) compo- 
nents condense simultaneously at T2 (solid line). For 23 Na 
atoms with antiferromagnetic interactions, double condensa- 
tions persist according to our theoretical study. The |+) com- 
ponent (denoted by +) condenses first, which is then followed 
by the condensation of the | — ) component (denoted by x). 
The |0) component is unpopulated in the low temperature 
limit. For 87 Rb atoms with ferromagnetic interactions, our 
study reveals the potential for triple condensations. First, 
the |+) component condenses (denoted by o), which is then 
followed by the second condensation for the | — ) component 
(denoted by *), and finally the third condensation for the |0) 
component occurs (denoted by o). 
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with the quantization axis taken along the z-axis direc- 
tion. It is easy to check that both the total number of 
atoms and the total magnetization 



N 
M 



dr(\i> + \ 2 + \^ \ 2 + \^\ 2 ), 
dr(h/>+| 2 -h/M 2 )= / dfVteylfc 



commute with the above Hamiltonian JBJ, and are thus 
constants of motion. To study the minimal energy 
ground state, we therefore introduce two Lagrange mul- 
tipliers /i and 77, to fix the total atom number and mag- 
netization of the system in our numerical minimization. 
It turns out that /_t is in fact the chemical potential of the 
system and 77 is an effective magnetic field. The Gibbs 
free energy is then given by 



G 



where 



H — fiN — rjM 

fWWy^^wV'i}, (7) 
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2m 



A 4 + Kxt(r) 



The atomic interactions are conveniently parametrized 
through the two s-wave scattering lengths eto and a-i be- 
tween two spin-1 atoms 0. [J 0] 
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In this study, we attempt to find the mean field ground 
state of our system, which corresponds to the state with 
the lowest Gibbs free energy. 



B. Hartree-Fock- Popov theory and the two-fluid 
model approximations 

The field operator ip(r, t) evolves in the Heisenberg pic- 
ture according to 

ih~r/>(?,t) = [i/>,G\. 

For the system of a spin-1 Bose gas as considered here, 
the above equation becomes 



3 

ifi^i>o(^t) = £00^0 + c ^2 (fjipj) Vt) + C2 (i>+i>+ + ip-ip-^j i>o + 2iplip+ip- , 

j 

ih^tjj_(f,t) = C — ^_ + 77^_ + cp ip- + c 2 + V'oV'o - + i>+foi>o 



(8) 



Following the standard mean field theory procedure, i.e. taking ip — <p + Sip with <p = (-0), after tedious manipulations 
and calculations, we obtain a set of coupled Gross-Pitaevskii equations for the superfluid components including their 
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interactions with the noncondensed atoms as 

h 2 



d 



-— V + Vext - M - V + c o(n + n + ) + c 2 (n + + n Q - n_ + n + ) 
2m 



t, 9 A 



-— V 2 + V ext - M + c o( n + n o) + c 2(n+ + n. 
Zm 
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ih 
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V 2 + V ext — (i + r) + c (n + n T ) + c 2 (n_ + n Q - n + + n^) 



and equations for Sipi, 

a 



ih-Q t H+{r,t) 



ih—5ip (r,t) 
d 
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2m 



V 2 + Kxt - A* _ V + c oi n + n+) + c 2 (2n + + n - ra_) 



— V 2 + Vcxt - M+ °o( n + no) + c 2 (n + + n_) 
Zm 



-— V 2 + F ex t - A* + V + co(n + n_) + c 2 (2n_ + ?i - n+) 
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where n = J2i n i = Si(l^i| 2 + n J) 1S the total den- 
sity of the atomic gas, with nf = (SipjSipi) the normal 
(noncondensed) gas density of the ith component. In- 
stead of the Hartree-Fock-Bogoliubov (HF B) approxima- 
tion as employed by Isoshima et al. jl9|, we have used 
the Hartree-Fock- Popov (HFP) approximation to obtain 
the above equations. Within the HFP approximation, we 
neglect terms proportional to the anomalous nonconden- 
sate density (SipiStpj) as well as their complex conjugates. 

We have also neglected terms proportional to (Sip} Sip j) 
for i 7^ j, similar to the random phase approximation. A 
more detailed formal discussion of the HFP theory can 
be found in Refs. 0, |2j|, and for the calculation of 
the phase diagram of Bose-Einstein condensation, it is 
an excellent approximation as confirmed recently in a set 
of detailed comparisons with experiments |2l| . In addi- 
tion, as will become clear later, the HFP approximation is 
also efficient from the numerical point of view, especially 
near regions of temperatures close to (but below) the 
critical temperature. The HFB approximation, on the 
other hand, is more difficult to handle numerically |19| . 
Although more rigorous at very low temperatures, the 
HFB approximation is expected to agree with the more 
transparent HFP approximation at higher temperatures. 
In deriving the equations for Stpi, terms proportional to 
Sipj, 6ipj, and Sipj for j ^ i are also neglected. This is 
equivalent to the neglect of the "hole" component in the 
HFB approximation, and is thus expected to have has a 
minor effect except very close to the zero temperature. 



In the HFP approximation we adopt here, the normal 
fluid for noncondensed atoms is determined through the 
semi-classical approximation. We thus take —iKV — ► p, 
and approximate its distribution by the standard Bose- 



Einstein distribution in the phase space of {p, r}, 

dp 1 



n T (r) 

lV> J {2-Khf e^v-.r)lh B T _y 
with the HFP single particle energy spectrum, 



(11) 
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~z h V ext - /i + c (n + n ) + c 2 (n + + n_) . 

Zm 



+ Vext - A* - V + c o(n + n+) 
n - n_) , 



P 

2m 

+c 2 (2n_ 



Kxt - A* + V + c o(n + ri-) 
n - n+) , 



(12) 



which are obtained by substituting Sipi (r , t) = 
exp[~iei(p,f)t/fi]ui(r) into Eqs. (fTU|) with it;(r) the 
eigenfunction for the excitation of the ith component. 

Thus, we have formulated a coupled set of equations 
for both the superfluid and the normal fluid; they are 
Eq. |J3J for the condensed part and Eqs. l(TT|) and ifT^ 
for the noncondensed atoms. These are the basis for our 
numerical investigations to be presented below. 



IV. NUMERICAL METHOD 

In our numerical studies, we follow a standard proce- 
dure and the following algorithm for the self-consistent 
solution of the coupled equations 10, (|11|) . and l|12|l as 
an extension of the single component gas studied earlier 
[2^|. Our algorithm is divided into the following steps: 

• We find the condensate wave function 4>i{r) and 
the chemical potential /i for a set of fixed normal 
gas density nf(r), by propagating Eqs. @ in the 
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• We compute the updated energy spectrum and nor- 
mal gas density nf(r) from Eqs. Hll|) and H12|) 
using the new condensate wave function and the 
chemical potential. 



• We normalize the total number of atoms to N and 
adjust r\ appropriately [26| . 

• We repeat the above steps until final convergence is 
reached. The convergence criterion is set to be that 
the condensate fraction N c (T) /N and the magne- 
tization fraction M/N of successive iterations differ 
by less than 10 -11 for most temperatures and less 
than 10 -5 near the phase transition temperature 
region. 

At temperatures higher than the first BEC transition 
point, the above procedure converges rather quickly as 
the superfluid component <f>i is none existent. In this case, 
we only need to solve Eqs. I|f f(l and (|f 2(1 sclf-consistently 
by adjusting fi and 77. 

V. RESULTS AND DISCUSSION 

In this study, we focus on the illustration of our theory 
for atoms inside a spherically symmetric harmonic trap 

V cxt (r) = ^niLU 2 r 2 . (13) 

We take N = 10 6 and w = (2tt)I00 Hz, and use the spin- 
1 atom parameters for 23 Na and 87 Rb atoms as given in 
Table HJ Clearly it is antiferromagnetic (02 > 0) for 23 Na 
atoms and ferromagnetic [c^ < 0) for 87 Rb atoms. 

TABLE I: Atomic parameters for 23 Na and 87 Rb atoms |2& 
[2S|. ao and a 2 are in units of Bohr radius and Co and C2 in 
units of 10" 12 Hz cm 3 . 





ao 


CL2 


CO 


C2 


23 Na 


50.0 


55.0 


15.587 


0.4871 


87 Rb 


101.8 


100.4 


7.793 


-0.0361 



A. Na atoms with antiferromagnetic interactions 

The interaction between 23 Na atoms is antiferromag- 
netic, i.e. C2 > 0. The phase diagram we obtain is shown 
in Fig. It clearly reveals the double phase transitions: 
one for the |+) component and the other for the |— ) 
component. The |0) component of the condensate never 
shows up because the antiferromagnetic interaction fa- 
vors an antiparallel alignment of the atomic spin, which is 
equivalent to a coherent superposition of the |+) and |— ) 
states as explained in the discussion of order parameter 




FIG. 2: Double condensations for a spin-1 gas of Na atoms. 
The upper left panel shows the total condensed fraction vs 
temperature and total magnetization. Similarly, the upper 
right one shows the fraction of condensed |+) component, the 
lower left the condensed |0) component, and the lower right 
the condensed | — ) component. 



symmetry at zero temperature in Ref. |3j . Our mean field 
result is also consistent with that of Isoshima et al. |19| . 
Similar to the case of an ideal gas, the transition tem- 
perature of the |+) component increases monotonically 
with M/N while that of the |— ) component monoton- 
ically decreases. When the temperature decreases, the 
first condensed component is |+) because M > 0; the 
second condensed component is |— ), which condenses at 
temperatures when N? + NT — NT > M. Figure shows 
typical density distributions of different components for 
a 23 Na gas. We see that |0+| 2 and |^-| 2 are always misci- 
ble |3Ctl3l| and distributed mostly near the central region 
of the trap. We also note that \<j)o\ 2 is always zero within 
this mean field study. All three components of the nor- 
mal gas coexist. Both and r\F_ peak at the edge of the 
condensate because of the shape of the net interaction 
potentials between the condensate and the normal gas. 
Uq is much flatter since \(j>o\ 2 is zero. 

We now comment on a particular feature related to 
the asymptotic behavior of the spin-1 gas of 23 Na atoms 
as T — ► for M = 0. The full quantum theory pre- 
dicts a ground as a superfragmented Fock state with 
atoms equally distributed among the three spin compo- 
nents \N+ = N/3, N = A/3, A_ = N/3) HHHH^. 
Such a state would give rise to a number fluctuation of 
order of A 2 , and is impossible within the present mean 
field treatment. The mean field ground state is known 
to be \N+ = N/2, N = 0, A_ = N/2) jHIH, consis- 
tent with our results. In an actual experiment, it is most 
likely that the mean field ground state is observed be- 
cause the full quantum state is not stable against various 
external sources of fluctuation or noises, e.g. that of the 
unshielded magnetic field |3^|, a small deviation of the 
total magnetization M from zero |32|. or a temperature 
being not exactly zero. The mean field ground state, on 
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FIG. 3: Typical density distributions for different spin com- 
ponents of a 23 Na gas (T/T c = 0.43, M/N = 0.4). The upper 
panel is for the condensate and the lower one for the noncon- 
densed atoms. 



the other hand, is more robust against these noise. In 
our numerical calculations, it is really impractical to set 
the temperature microscopically close to but above zero 
to probe the real ground state (for M = 0). We therefore 
enforced the ground state structure such that it asymp- 
totically approaches that of the mean field ground state 
with a decrease of the temperature as shown in Fig. [21 
With this convention, an related issue arises: the equiv- 
alence between states \N+ = N/2, N = 0, 7V_ = N/2) 
and \N+ — 0, No = N,N^ = 0) at zero temperature as 
first pointed out by Ho We note, however, that this 
equivalence is based on the assumption of an environ- 
ment perfectly free of magnetic fields. The presence of 
even a tiny magnetic field, which is inevitable in the real 
world, will destroy this equivalence and causing the real 
ground state to be \N+ = N/2,N = 0, = N/2), a 
convention we chose as indicated in Fig. [21 



B. Rb atoms with ferromagnetic interactions 

The phase diagrams for 87 Rb atoms with ferromagnetic 
interactions (02 < 0) are shown in Figs. [Hand El Only a 
sparse set of points weas made available in the early work 
of Isoshima et al. because the numerical solution be- 
comes far more difficult to converge in this case. Based 
on our results, we see that when the temperature of the 
system decreases, triple condensations occur in general. 
When M > 0, the first condensed component is the |+) 
state (for T < Ti), the second one is the |— ) state (for 
T < T 2 ), and the last one is the |0) state (for T < T 3 ). 
Our results show that the |+) component first condenses 
at T\ and its population increases with decreasing tem- 
perature until T2, at which the |+) condensed component 
is a little more than the total magnetization M. When 




FIG. 4: The same as in Fig. [5] but for a gas of 87 Rb atoms. 



temperature is lower than T2, the |— ) component begins 
to condense as well. The two condensed components in 
states |+) and |— ) both increase with decreasing temper- 
ature until the third critical temperature T3, at which 
the |0) component starts to condense. Once the |0) com- 
ponent condenses, the |— ) component starts to decrease 
and becomes very close to zero, while the |+) component 
is almost constant with respect to further decreasing of 
the temperature. This trend continues until the tempera- 
ture is lower than T4, when the condensed |0) component 
starts to decrease with decreasing temperature while the 
populations of the |+) and |— ) condensed components in- 
crease. For the special case of M = 0, on the other hand, 
we observe only double condensations; the |0) component 
condenses first, followed by the simultaneous condensa- 
tion of both the |+) and |— ) components. This is again 
due to the special symmetry requirement that the |+) 
component must be the same as the |— ) component in 
order to keep M = 0. We note that in this case the 
fraction of condensed |0) component can reach as high as 
94% at finite temperatures, much higher than the ~ 50% 
at zero temperature. 

Figure[fJldisplays typical density distributions for a gas 
of 87 Rb atoms at different temperatures for M/N = 0.6. 
The right column corresponds to T € (T3, T2], where only 
the |+) and |— ) components are condensed and the |— ) 
component is quite small and is spatially located at the 
edge of the |+) component. Quite generally, we note that 
with the condensation of a component, its corresponding 
normal gas component would have a lower density. For 
instance, the normal gas density of the |+) component is 
low in the trap center where the |+) condensed compo- 
nent resides. The middle column of Fig. [SJis the typical 
density distribution when T € (T±,T^. The condensed 
|+) component stays at the center and is surrounded by 
the |0) component. The condensed |— ) component is too 
small to be visible directly, but can be perceived from the 
shallow well in its normal gas component, which indicates 
that the condensed |— ) component is not zero and is lo- 
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FIG. 7: The lowest excitation level for a gas of 87 Rb atoms at 
M/N = 0.6, T/T c = 0.34 (right before the condensation of the 
|— ) component). The inset shows the details of a zoomed-in 
plot near the minimum. 



FIG. 5: Double condensations for Rb atoms when M = 
(the upper panel) and triple condensations for M/N = 0.6 
(the lower panel). The solid line denotes the fractional pop- 
ulation of the condensed |+) component, the dot-dashed line 
denotes the jO) component, and the dashed line denotes the 
| — ) component. 
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FIG. 6: Typical density distributions for a gas of 87 Rb atoms 
at M/N = 0.6. The left column is T/T c = 0.11, the middle 
one is 0.21, and the right one 0.32. The notations are the 
same as in Fig. El 

cated around the edge of the condensed |+) component. 
The left column of Fig. [fllshows the density distributions 
when T 6 (0, T4], where all three condensed components 
coexist near the center of the trap and are surrounded by 
their normal gas components. 

These results for 87 Rb atoms can be understood in 
terms of the interplay of three factors: ferromagnetic 
atom-atom interaction (02 < 0), the M conservation, and 
the miscibility between and among different components. 
The ferromagnetic interaction favors the most populated 
state, the M conservation sets an upper limit on the frac- 



tion of the condensed |+) component, and the immiscibil- 
ity between the condensed |+) and |— ) component sets an 
upper limit on the total fraction of the condensed |+) and 
|— ) components. For instance, in the region T 6 {T 2 , Ti], 
only the |+) component condenses. The ferromagnetic 
interaction plays a dominant role and thus more atoms 
condense into the |+) state with decreasing temperature. 
In the region of T G (T3 , T2] , the M conservation and the 
immiscibility begin to take their effect. The M conser- 
vation causes the increases to the condensed |+) and |— ) 
components to be almost identical, while the immiscibil- 
ity makes the condensed |— ) component stay outside the 
condensed |+) component. The system becomes unstable 
with the increase of the |+) and |— ) components because 
with more condensed |+) component, the condensed |— ) 
component must be pushed out further. Near the third 
critical temperature T 3 , the condensed |— ) component 
suddenly decreases to almost zero and the condensed |+) 
component decreases to about M. Approximately, the 
total decreased amount from the |+) and the |— ) com- 
ponents becomes the condensed |0) component. The sys- 
tem enters the region T £ (Ti, T3] in which the condensed 
|0) component increases steadily with lowering tempera- 
ture because the condensed |+) and |0) components are 
miscible. The condensed |+) component is almost in- 
dependent of the temperature because of the M conser- 
vation. With decreasing temperatures, more and more 
condensed |0) component finally suppresses the immis- 
cibility between |+) and |— ) component at T4, all three 
condensed components become miscible, and both the 
|+) and |— ) components increase to keep M conserved 
while the |0) component decreases. 

Figure shows the lowest excitation energy e = e(p = 
0) for the three components of a 87 Rb gas at M/N = 0.6, 
T/T c = 0.34 (right before the condensation of the |— ) 
component). We see that the energy for the |— ) compo- 
nent is lower than the corresponding ones for the other 
two states and takes a minimum near the spatial location 
of r = 9 /im, which is at the edge of the condensed |+) 
component. This result confirms that the |— ) compo- 
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nent condenses before the |0) component and surrounds 
the condensed |+) component. 

VI. CONCLUSION 

We have studied the thermodynamics of Bose-Einstein 
condensation for a gas of spin-1 atoms with ferromag- 
netic and antiferromagnetic interactions using the mean 
field Hartree-Fock-Popov theory and the semiclassical ap- 
proximation for the noncondensed components. Our re- 
sults show that for antiferromagnetic interactions, dou- 
ble phase transitions persist as in a noninteracting gas: 
when M > 0, first the |+) component condenses, which 
is followed by the condensation of the |— ) component on 
further decreasing of the temperature. The |0) compo- 
nent never condenses. For ferromagnetic interactions, on 
the other hand, our calculations reveal that the phase 
diagram becomes more complicated and a triple conden- 
sation scenario arises with decreasing temperatures: first 
the |+) component condenses, which is followed by the 
second condensation of the |— ) component, and the third 
one for the |0) component. When the |+) and |— ) compo- 
nents are the only condensed ones, they are immiscible. 
When all three components condense and the tempera- 
ture is lower than T4, they become miscible because of the 
presence of a large condensed |0) component. We have 
compared the numerically computed transition tempera- 
tures with that of an ideal gas as in Fig. ^ An overall 



lowering of the various transition temperatures due to 
atom-atom interactions is seen, consistent with the case 
of a single component interacting Bose gas, where the 
interaction-caused shift to the transition temperature has 
been actively studied [35| ■ Quite generally a repulsive in- 
teraction tends to lower the transition temperature for a 
single component Bose gas [24j. In the case of a spin-1 
Bose gas considered here, Co > and c ^> |c 2 | consti- 
tutes an overall repulsive interaction. 

Finally, we note there also exists the possibility of a 
ferromagnetic phase transition for 87 Rb atoms, in addi- 
tion to the Bose-Einstein condensation as studied here. 
In fact, as was investigated recently by Gu and Klemm 
|3fil |. the ferromagnetic transition is generally predicted 
to occur before, i.e. at temperatures higher than, the 
Bose Einstein condensation. The present study, how- 
ever, remains unchanged because we treated the system 
within the global constraint of the conservation of total 
magnetization, distinct from that required for a sepa- 
rate ferromagnetic phase transition [36j . As is evidenced 
from recent experiments, the total magnetization M is 
well conserved, even better than the conservation of the 
total number TV [Tlj . 
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